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We consider various approaches to treat the phases of a qutrit. Although it is possible to represent 
qutrits in a convenient geometrical manner by resorting to a generalization of the Poincare sphere, we 
argue that the appropriate way of dealing with this problem is through phase operators associated 
with the algebra su(3). The rather unusual properties of these phases are caused by the small 
dimension of the system and are explored in detail. We also examine the positive operator-valued 
measures that can describe the qutrit phase properties. 
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I. INTRODUCTION 

The emerging field of quantum information, which em- 
braces areas of futuristic technology such as quantum 
computing, quantum cryptography, and quantum com- 
munications, has been built on the concepts of entan- 
glement and qubits P, . The full appreciation of the 
complex quantal properties of these two ideas has pro- 
vided powerful physical resources for new schemes that 
herald results that cannot be achieved classically 3]. 

Recently, the exploration of higher dimensional quan- 
tum systems has finally received the attention it rightly 
deserves. One could think that this represents a mere 
digression in a hot topic. However, qutrits have several 
interesting properties worth exploring the efficiency 
and security of many quantum information protocols are 
improved using qutrits ^^ii £i 1] , and larger violations of 
nonlocality via Bell tests are expected to occur for sys- 
tems of entangled qutrits • 

In the modern parlance of quantum information the 
concept of phase for a qubit (or a qutrit) is ubiquitous. 
However, this notion is rather imprecise. Phases for 
three-level systems have been handled by invoking fuzzy 
concepts such as the phase of the associated wavefunc- 
tion Sometimes, the problem is reduced to the opti- 
mal estimation of the value of the phase shift undergone 
by the qutrit 'Tl'|. 

When comparing phases of two states, it is usually 
assumed that the relative phase is obtained from the 
argument of their inner product. In this perspective, 
the phase is considered as a state parameter. In recent 
years, we have learned that this relative phase shift can 
be of various origins, namely, it can be purely dynam- 



ical or purely geometrical or both. Presently, there is 
an immense interest in geometric phases in quantum op- 
tics [T^ Il3l[ , especially in connection with quantum com- 
puting applications "^l- In fact, these phases are linked 
to the geometry of the state space: for a qubit, this space 
is the coset space SU(2)/U(1), the well-known Poincare 
sphere, while for a qutrit, a geometrical picture of the 
corresponding generalization to SU(3)/U(2) has been re- 
cently presented [l5|. 



We emphasize that these notions, though well estab- 
lished in the classical limit, are not easily extrapolated 
into the realm of the quantum world. Since the phase 
is a physical property, it must, in the orthodox picture 
of quantum mechanics, be associated with a selfadjoint 
operator or at least with a family of positive operator- 
valued measures (POVMs). In this spirit, phase opera- 
tors for the algebra su(2), which describes qubits, have 
been previously worked out [l^ ll^ ITsL IT^ . as well as the 
optimal POVM for this problem [23]. The main goal of 
this paper is to work out a nontrivial extension to su(3) 
of the results available for su(2), enabhng us to introduce 
phase operators for qutrits with a clear physical picture. 
This seems of such fundamental importance, that it is 
surprising that such a task has not been undertaken long 
time ago. We thus trust that this will be of relevance 
to workers in the various experimental fields currently 
under consideration for quantum computing technology 
and in quantum optics, in general. 
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II. POINCARE SPHERE FOR A QUTRIT 

We first briefly recall the salient features of the 
Poincare sphere representation for a qubit, with a view 
of preparing its generalization for a qutrit along the lines 
of Ref . j3| • A. qubit lives in a two-dimensional complex 
Hilbert space 7i^^^ spanned by two states: |1) and |2). To 
get a useful parametrization of the state space of a gen- 
eral qubit described by the density matrix p, we observe 
that the Pauli matrices da together with the identity 1 
form a complete set of linearly independent observables, 
and that any selfadjoint (trace class) operator can then 
be written as 

p=i(i + n.^). (2.1) 

The physical condition /5 > holds only when |n| < 1. 
Hence, the state space coincides with the Bloch ball, and 
the set of pure states (p^ = p) with the boundary of this 
ball |n| = 1, which is the Bloch sphere . The general 
pure state 

I*) = sin(6l/2)|l) + e"^ cos(6l/2)|2), (2.2) 

is represented by the unit vector 

n= (sin ^? cos 0, sin0 sin (/), cos 0). (2-3) 

The angle d is obviously related to the proportion of |1) 
and 1 2) in the composition of the state, while the param- 
eter is routinely interpreted as the quantum phase as- 
sociated with the qubit and canonically conjugate to the 
inversion 0|. We note in passing that diametrically 
opposite points on correspond to mutually orthogonal 
vectors in liS^'> . 

Consider now a qutrit, living in a three-dimensional 
complex Hilbert space Ti*^'^^ spanned by |1), |2), and |3). 
The roles of SU(2) and the Pauli matrices are now played 
by the group SU(3) and the eight generators of the cor- 
responding su(3) algebra. A convenient set of Hermi- 
tian generators are the Gell-Mann matrices [23, A,- 
(r = 1, . . . , 8), which obey the commutation relations 

[A„A,] = 2i/,,tAt, (2.4) 

where, above and in the following, the summation over 
repeated indices applies. The structure constants jrst 
are elements of a completely antisymmetric tensor spelled 
out explicitly in Ref. |15j|, for example. 

A particular feature of the generators of SU(3) in the 
defining 3x3 matrix representation is closure under an- 
ticommutation 2^ 

- , 4 
{A^, AJ = -(5,.,1 + 2drst\t, (2.5) 

where now drst form a totally symmetric tensor. 

For the following, a vector-type notation is useful, 
based on the structure constants. The / and d sym- 



bols allow us to define both antisymmetric and symmet- 
ric products by 

(AAB), = j,,tA,Bt = -{^KK\ 

(2.6) 

(A*B)r = VSdrstAsBf = -|-(B* A)r. 

Given a density matrix p we can expand it in terms of 
the unit matrix 1 and the A,- in the form 

/5-i(l + V3n.A). (2.7) 

This is the equivalent to the Bloch ball for a qutrit. For 
a pure state the analogous Bloch sphere is defined by the 
condition 

n-n = l, n*n = n. (2-8) 

Thus, each pure qutrit state corresponds to a unique unit 
vector n g 5^, the seven-dimensional unit sphere. In 
addition, this vector must obey the condition n ★ n = n, 
which places three additional constraints, thus reducing 
the number of real parameters required to specify a pure 
state from the seven parameters needed to specify an 
arbitrary eight-dimensional vector to four. 

In view of our discussion for qubits, it is clear that 
normalization and a choice for the arbitrary overall phase 
allow us to write these four parameters for any pure state 
as 

l^-) = sin(^/2)cos(e'/2)|l) 

e*=sin(C/2)sin(6l/2)|2) 
+ e^^"cos(^/2)|3), (2.9) 

Again, Q and f determine the magnitudes of the compo- 
nents of 1 4'), while we can interpret (/)i2 as the phase of 
|1) relative to |2) and analogously for 0i3. We can easily 
obtain the expressions for n in these local coordinates. 

Some interesting geometric properties of this Poincare 
sphere are discussed in Ref. 4J. In particular, it is easily 
seen that for two unit vectors n and n' representing pure 
states 

27r 

< arccos(n • n') < — , (2.10) 
3 

so mutually orthogonal vectors in H!^^^ do not lead to an- 
tipodal or diametrically opposite points on the Poincare 
sphere, but to points with a maximum opening angle of 
27r/3. 

III. PHASE OPERATORS FOR A QUTRIT 

Although there is a widespread usage of dealing with 
the qutrit phases as state parameters, this is not an or- 
thodox way of proceeding, according to the very basic 
principles of quantum mechanics. 
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To gain further insights into this obvious although al- 
most unnoticed point, we stress that the complete de- 
scription of a qutrit involves the nine operators 



5,, 



(3.1) 



where |i) is a basis vector in Ti'^^. The three "diagonal" 
operators Su measure level populations, while the "off- 
diagonal" ones Sij represent transitions from j to level i. 
One can easily check that they satisfy 



[Sij, Ski] — SjkSii — SiiS^ 



(3.2) 



which are the commutation relations of the algebra 
u(3) 24]. 

Because of the trivial constraint 5ii + S22 + S33 = 1, 
only two populations can vary independently. For this 
reason, we shall work with two independent traceless op- 
erators 



^12 - o('5'22 - S'll), S23 - - 



(^33-^22), (3.3) 



that measure atomic inversions between the correspond- 
ing levels. In atomic systems, the selection rules usually 
rule out one of the transitions and therefore the two inde- 
pendent inversions are automatically fixed. For a general 
qutrit, these inversions can be arbitrarily chosen. 

The commuting operators Sf2 and 6*13 constitute a 
maximal abelian subalgebra for the qutrit (known as Car- 
tan subalgebra). From the discussion of the previous 
section, we expect 8^2 and S23 to be conjugate to the 
corresponding (independent) phases of the qutrit. 

Note that {812,8^2) and (S'23,5'23) correspond to the 
qubits 1 4-> 2 and 2 <-> 3. However, these two qubits are 
not independent, since Eq. (|3.2|l imposes highly nontriv- 
ial coupling between them. 

At the operator level, the equivalent to the decomposi- 
tion of a complex number in terms of modulus and phase 
is a polar decomposition [25j. S ince tS'21 — '^121 seems 
appropriate to define 



'S'12 — R12 Ei2j 



(3.4) 



where the "modulus" is R12 = ^8 12 5*21 and E12 — 
exp(i0i2), (/'12 being the Hermitian operator represent- 
ing the phase. 

One can easily work out that a unitary solution of 
Eq. H3.4|l is given, up to an overall phase, by 



^i2 = |l)(2|+e^^" |2)(1| 



|3)(3|, 



(3.5) 



where the undefined factor e"^° appears due to the uni- 
tarity requirement of E12. The main features of this op- 
erator are largely independent of (/jq, but for the sake of 
concreteness, we can make a definite choice. For exam- 
ple [T9|. for a qubit defined by a linear superposition of 
the states |1) and |2), the complex conjugation of the 
wave function should reverse the sign of 0i2, which im- 
mediately leads to the condition e"^" = — 1. We conclude 



then that a unitary phase operator that preserves the 
polar decomposition of Eq. 13.411 can be represented as 



E12 



|1>(2|-|2)(1|-|-|3)(3|. 



(3.6) 



The eigenstates of ^12 are those of Ei2, and easily found 
to be 

l'/'?2> = |3), |</.±2> = ^(|2)±*|1)), (3.7) 

with the corresponding eigenvalues of (f)i2 and ^1^/2, 
respectively. This is a remarkable result. It shows that 
the eigenvectors \(l)i2) ^ook like the standard ones for a 
qubit. However, the "spectator" level |3) is an eigenstate 
of this operator, which introduces drastic changes. In 
other words, the phase of the qubit 1 2 "feels" the 
state 1 3). 

An analogous reasoning for the transition 2^3 gives 
the corresponding operator i?23 



^23-|2)(3|-|3)(2| + |l)(l| 
with eigenvectors 



ID 



10: 



^(|3)±*|2)), 



(3.8) 



(3.9) 



and the same spectrum as before. 

As for the operator i?i3, one must be careful, because 
it connects the lowest to the highest vector. In fact, the 
polar decomposition in this case gives as a unitary solu- 
tion 

^3 = a|3)(2|-6*|3)(l|+6|2)(2|+a*|2)(l| + |l)(3|, (3.10) 



with the condition 



|6p — 1. There are also nonuni- 



tary solutions to the polar decomposition, but they lack 
of interest to describe a phase observable in our context. 

Note that the general solution (|3.1U|) has the desirable 
property -Ei3|3) = |1). On physical grounds, we argue 
that the state |2) should be a "spectator" for the transi- 
tion 1 3. Thus we impose i?i3|2) oc |2), which is only 
possible if a = and we have that 



i?i3 = |l)(3|-|3)(l| + |2)(2|, 
with eigenvectors 

|2), |0±3) = ^(|3)±*|1)). 
With this choice we are led to 



^?3) 



E12E23 El3, 



(3.11) 



(3.12) 



(3.13) 



which clearly displays the quantum nature of this 
phase ,27]. Note, in passing, that 



[i?12,i?23] = [£'23,^12] = 0, 



(3.14) 



and [R23,Ri2] — 0, so the interference between differ- 
ent channels (i.e., the noncommutativity of S'12 and 823) 
is due to the nonncommutativity of the corresponding 
phases. 
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IV. POSITIVE OPERATOR MEASURES FOR 
THE QUTRIT PHASES 

The unusual behavior exhibited by the description of 
qutrit phases in terms of Hermitian operators can be 
considered to some extent exotic. One may think it 
preferable to represent qutrit phases by using a positive 
operator-valued measure (POVM) taking continuous val- 
ues in a 2n interval. 

We briefly recall that a POVM [2^ associated to an ob- 
servable (/) is a set of linear operators A(0) (0 < (/) < 27r), 
depending on the continuous parameter (p and furnish- 
ing the correct probabilities in any measurement process 
through the fundamental postulate that 



P((/.) = Tr[/5 A(0)]. 



(4.1) 



The real valuedness, positivity, and normalization of 
P{(j)) impose 



At(0) = A(0), A(</.)>0, 



2-K 



d(j) A((/)) = i. (4.2) 



where the integral extends over any 27r interval of the 
form ((/iQ, 00 + 27r), (pQ being a fiducial or reference phase. 
Note that, in general, A(0) are not orthogonal projectors 
as in the standard von Neumann measurements described 
by selfadjoint operators. 

From our previous discussion, it is clear that we ex- 
pect some c omp lementarity between phases and inver- 
sions ^3, 130, . If we observe that 



-'^'/2|l)(l|+e^'^72|2)(2| + |3>(3|, (4.3) 



and argue that phase-shift operators must be 27r peri- 
odic, we impose that any POVM A(0i2, ^23) for a qutrit 
should satisfy 



e'20'S^.A(0,2,(/.23)e-*''^'^- 



+ </''>23), 



(4.4) 



A( 



2, 023 + 0")- 



One can work out that the general POVM fulfilling these 
requirements must be of the form 



A( 



1 



[7126 



i(2(/)i2— 023) 



|2)(1| 



^23e'(2te-'^i2) (2| + 7i3e*('^i=+'^^^) |3) (1| 
h. c.]}, (4.5) 



where h. c. denotes Hermitian conjugate, jij < 1 are 
real numbers and (jjij is the relative phase between states 
\i) and \j). These relative phases coincide precisely with 
the polar part of the realization of su(3) on the torus 
constructed in Ref. [23|. If we chose the 7^^ different, 
say 712 — 1 and the other two below the unity, then the 
expectation value of this POVM could reach the value 
zero for the superposition states (|1) -I- exp(i6')|2))/V2. 



However, for superpositions of states |1) and |3) or |2) and 
1 3), the expectation values of the POVM would always 
be greater than zero. Since there is no physical reason to 
assign special relevance to one specific superposition of 
the states, we assume that the POVM must be symmetric 
with respect to the states, which leads to 



7 = 712 = 723 = 713- 



(4.6) 



Moreover, we make henceforth the choice 7=1 because 
only for this choice the POVM can attain the expectation 
value zero for some particular state. 

In contrast with the result of Eq. (|3.13|) formulated 
in terms of operators, now there are only two relevant 
phases in the qutrit description: the third can be inferred 
from the other two, as in the classical description. 

The proposed POVM provides qutrit phases where any 
values of 0i2 and 023 sue allowed. However, note that the 
probability density induced by this POVM can be written 
as 

^2,023) = J^{1 + [pi2e'^^'^''-'^'''> + p23e'^^'^''-'^'''> 
„«(012 + fe) c.]}. 



Pi3e 



(4.7) 



where pij = {i\p\j) and c. c. denotes complex conjugate. 
Therefore, this continuous range of variation is not effec- 
tive in the sense that the values of P(0i2 , 023) at every 
point (012,023) cannot be independent, and we can find 
relations between them irrespective of the qutrit state. 
In other words, the complex parameters Cy can be deter- 
mined by the values of /'(0i2 , 023) at six points. Discrete- 
ness is inevitably at the heart of the qutrit phase . 

Finally, we shall consider a remarkable example of 
POVM particularly suited to describe the qutrit phase. 
We recall that for a single-mode quantum field, a POVM 
for the field phase can be defined in terms of radial inte- 
gration of quasiprobability distributions obtained using 
a coherent-state representation, much in the spirit of the 
classical conception js^, 13 • The natural generalization 
of this procedure to the qutrit problem involves the use 
of su(3) coherent states. In the Appendix we summarize 
the essential ingredients needed for this paper. Coherent 
states of a single qutrit are of the form 



1 



= (|3>+a|2)+a/3|l)). 



(4.8) 



where a and (5 are complex numbers and the normaliza- 
tion constant is 



Ca/3 - 1 + |a|'(l + 



(4.9) 



These coherent states generate a POVM over the qutrit 
state space via the projectors ja, /3) (a, /3|. 

As shown in the Appendix, the phases of a and /3 are 
just those of (a, /3| 5*3210, /3) and (a, /3|5'2i |a, /3), respec- 
tively, while the phase associated to (a, (i\S'i\ |a, /?) is just 
the product of the other two. Let us write 



a 



(3 = ri2e* 



(4.10) 
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and integrate the projectors |a, (3) {a, /3| radially over ri2 
and r23, with respect the measure [see Eq. (|A.11|) ] 



[l + |aP(l + |/3P)]3 
After some calculations one obtains 



(4.11) 



A(0i2,023) = ^{i + ^[e^(2^--*-)|2)(l| 



+ e 



i(2023-<;''12) 



|3)(2| 



+ e'('^^^+*^^)|3)(l|+h.c.]}, 



(4.12) 



which is just a specilized form of Eq. ()4.5|) and whose 
physical meaning is now clear. 



V. CONCLUDING REMARKS 

In this paper we have looked for possible descriptions 
of qutrit phases. Although it is possible to construct an 
extension of the Poincare sphere to qutrits, the orthodox 
way of dealing with any observable is to represent them 
by self adjoint operators. In this spirit, we have investi- 
gated a description of qutrit phases in terms of a proper 
polar decomposition of its amplitudes. Perhaps the most 
striking consequence of this description is that phases are 
discrete and do not commute. 

We have also considered alternative generalized de- 
scriptions in terms of POVMs. In these descriptions, 
phases appear as parameters rather than operators. Ad- 
ditivity of phases follows from the commutativity of the 
Cartan elements. Although these POVMs reflect some 
desirable properties of the classical phase, they show an 
effective discreteness, even if in principle a continuous 
range of variation is assumed. 
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APPENDIX: SU(3) COHERENT STATES 

In this Appendix, we briefly summarize the essen- 
tial ingredients of the construction of coherent states for 
three-level systems 35j. For concreteness, we shall con- 
sider fully symmetrical states of N three-level systems. 
In the Fock representation, we denote by |ni, n2, na) the 
state in which there are ni systems in level 1, n2 systems 
in level 2 and na systems in level 3. We observe that all 
these states can be generated from |0, 0, N) by repeated 



application of the usual collective operators S23 and 5*12 
[note that they coincide with H3.1|) . introduced for one 
qutrit, when iV = 1] as 



(^12)" (523)'"|0,0,iV) = 



Nlml 



{N - to)! y (to - n)\ 
X |n, TO — n,iV — to) (A.1) 

with < ri < TO < A^. Note that this is a simple ex- 
tension of the relevant formula for the two-level case. In 
analogy with the atomic coherent states for su(2), we 
define coherent states for qutrits as 

!«,/?) = v^e''^'' e"^'lO,0,A^), (A.2) 

where a and /3 are complex numbers and Maji is a nor- 
malization constant that we shall write as 



Q/3 



(A.3) 



where we have introduced the real quantity 

C^p^\ + \a\\\ + W). (A.4) 
In the Fock basis these states can be recast as 



|a,/3) = v^A^ ^ 

'i<n<m<N 

X a"^l3"'\n,m — n, N — m) 



\ 1/2 / \ 1/2 



(A.5) 



After some calculations one gets the following mean val- 
ues 



= {Q;(i\Szz\a,l3) 



N 

Cap ' 

N 



«2 = {a,P\S22\a,P) = (A.6) 

fil = {aP\S,^\a,P)^^\a\'\P\\ 

which immediately shows that the ratios of the average 
population numbers are given by 

na : n2 : rli = 1 : : \afm^. (A.7) 

On the other hand, one can also compute 

N 

(a,/3|S'32|a,/3) = 7;—") 

{a,(3\S2i\a,f3) = ^\a\^P, (A.8) 
N 

{a,(3\S3i\a,f3) = -p;—a(3. 

(--a/3 

The phases of a and /3 are then just those of 
(a, /3|S'32|a, /?) and (a,/3|S'2i|a,/3), respectively. Note, in 
passing, that the third phase associated to (a, PlS^i [a, /?) 
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is just the product of the other two, as it happens in clas- 
sical physics. 

The atomic coherent states with different amplitudes 
are not orthogonal 



(A.9) 



(ai,/?l|Q!2,/32) = 

[l + aW2(l + /?i/?2)]^ 

[1 + |ai|2(l + |/3iP)]^/2[l + |a2|2(l + |/32|2)]iv/2' 

but form an overcomplete set. In fact, it is easy to verify 
the following resolution of the identity 



{N + l){N + 2) 



j d/i \a,j3){a,f3\ = 1, (A.IO) 



where the measure dfi is 



[1+H2(1+|/3|2)]3 



(A.11) 



The above discussion pertains only to the fully symmet- 
ric subspace. Of course, it is enough for our purposes, 
although is can be extended to other subspaces in a very 
direct way. 
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